A quasi-exactly solvable model: two charges in a magnetic field, 
subject to a non-Coulomb mutual interaction 
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Abstract 

We extend the class of QM problems which permit for quasi-exact solutions. Specifically, we consider 
planar motion of two interacting charges in a constant uniform magnetic field. While Turbiner and 
Escobar-Ruiz (2013) addressed the case of the Coulomb interaction between the particles, we explore three 
other potentials. We do this by reducing the appropriate Hamiltonians to the second-order polynomials in 
the generators of the representation of SL(2, C) group in the differential form. This allows us to perform 
partial diagonalisation of the Hamiltonian, and to reduce the search for the first few energies and the 
corresponding wave functions to an algebraic procedure. 


^ Introduction 

o 

lO In quantum mechanics, quasi-exactly solvable (QES) problems are those which permit analytic derivation 
a hnite part of the (generally, inhnite) spectrum. They allow for partial reducing of the Hamiltonian 
a block-diagonal form. In these cases, one can find nonperturbative solutions for Hamiltonians of some 
Opecial types. These special types are obtained from Hamiltonians of broader (generally, unsolvable) classes 
J^y constraining the entering parameters. 

^ Exploration of these systems was pioneered by Turbiner [1] . Later, this kind of systems was independently 
• rediscovered by Kamran and Olver [2]. A different approach (often termed analytic) was developed by 
^^shveridze [3]. Various approaches to QES problems and some advances along this line of research are 
qf)resented in |1]. Among the recent advances in this field, notable are the discovery of ‘semiclassically quasi- 
exactly solvable potentials’ |S], as well as the results on studying elliptic and BCn models PE!. In the 
current paper, we employ a so-called algebraic approach, a detailed description whereof can be found in |8]. 

All exactly and quasi-exactly solvable one-dimensional Hamiltonians known so far can be reduced to 
second-order polynomials in the generators of a representation of the SL{2, C) group: 

h = ao^pj:j^ + bJ: , 

where J" are implemented as differential operators. Such Hamiltonians’ eigenvectors and eigenvalues are 
found through an algebraic procedure. Our goal is to broaden the class of problems which can be cast into 
the said form and which can, therefore, be treated by algebraic methods. A list of QES Hamiltonians of the 
form 

« = + nu (1) 
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was provided in [T] . When one considers the radial part of the eqnation for the relative motion in the problem 
of two particles in two or three dimensions, other QES Hamiltonians emerge. In distinction from the type (IT]) , 
they contain the hrst derivative. In [9] , several particular analytic solutions were described for a system of two 
electrons (interacting with Coulomb potentials) in a homogeneous magnetic field and an external oscillator 
potential. In [10], analytic solutions were found for two particles of opposite charges in a homogenous magnetic 
field. Both sets of solutions were obtained my means of a recurrence relation on the coefficients a„ in an anzats 
( Yin ) exp(...) for the wave function. In [TT], the algebraic structure of the Hamiltonian of a system of 
two electrons in an external oscillator potential was revealed. 

In this paper, we consider planar motion of two interacting charged spinless nonrelativistic particles subject 
to a constant uniform magnetic field. While Turbiner and Escobar-Ruiz [121 ESI E] addressed the case of a 
Coulomb interaction between the particles, we here explore other potentials that can permit for the afore- 
described reduction of the Hamiltonian. Of a special interest are the potentials resembling effective interaction 
between charges in realistic physical settings |14] . 

In Section |2l we briefly discuss the underlying idea of the method, an idea that ensues from the basic 
properties of representations of Lie groups. In Section [3l we provide the statement of the problem and carry 
out some prelusory calculations. In the subsequent sections we derive and explore three reducible potentials 
never presented in the literature hitherto. 


2 Basics of QES 


Let Jn be the generators of an irreducible finite-dimensional representation of the SL{2, C) group in some 
vector space: 

d„“.-hl = 4+ . l*Ll = -L . IL,41 = 2^“ , (2) 

where n = 2^' -|- 1 is the dimension of the space, and j is the spin of the representation. These relations are 
satisfied by the operators defined as 


Jn = ‘^jp - P dp , 
Jn ~ dp 5 

Jl = -j + pdp • 


(3) 


For each n, we find the corresponding vector space to be 


l/ = span{l,p, . (4) 

Indeed, these monomials are eigenvectors of J°, while and J~ are the raising and lowering operators. None 
of these operators can take us outside the representation space V. As a result of this, V stays invariant under 
the action of any polynomial function O of the operators Jn '■ 

dv<zv . (5) 


Solving the associated spectral problem 


ofip) = -^fip) 

over the vector space ([1]) is then equivalent to diagonalisation of an 
if one implements the basis vectors (0]) and the differential operators 
correspondingly. 
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( 6 ) 

n X n matrix. This becomes obvious 
(1^ with column vectors and matrices. 


In what follows, we shall be interested in operators of the type 

h = + bJ: , (7) 

because all exactly and quasi-exactly solvable one-dimensionsl Hamiltonians known hitherto can be reduced 
to this form. In this way, the spectral problem for the operator h , 

aapJnJn + fip) = “ /(p) , (8) 

is solvable by purely algebraic means. 


3 Formulation of the problem 

We consider a system of two non-relativistic charged spinless particles, (ei,mi) and ( 62 , 7712 ), moving on a 
plane orthogonal to a magnetic held B = Bz. We limit ourselves to a constant and uniform magnetic held, a 
restriction which permits us to use the symmetric gauge Ai ^2 = \B x pi^ 2 - The latter circumstance simplihes 
the problem greatly and allows us to write the Hamiltonian in the form of 

where h = c = 1, while pi 2 and pi 2 = —fV are the coordinates and momenta of the particles. 

It is natural to switch to the centre of mass reference frame (CMF): 

R = piPi + P2P2, P = Pi-P2 • (10) 

Canonically conjugated to them are the CMF and relative momenta: 


P = Pl+P2, P = P2P1-P1P2 , 


where 


M = mi + 1712 , 


M 


Using these formulae, one can prove that 


P = -iVi 


P = -fV, 


The above new coordinates and momenta enable us to dehne the CMF and relative angular momenta: 

L = R X P , i = p X p . (14) 

By direct calculation, one can check that the total pseudomomentum K and the total angular momentum 
are the integrals of the system: 

K = pi + eiAi+p2 + e2A2 = P+ qAR + ecAp , (15) 

= Pi X Pi + P 2 X P 2 = -T + , (16) 

[k,n] = [L^,n] = Q , (17) 


where 


_ , _ _ A e2 

g — 61 + 62 , Cc — P 2^1 — PlS -2 — TTlr- 

\mi m 2 
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ec being a coupling charge and nir being the reduced mass of the system: 


rur = 


mim2 


M 


After a unitary transformation 
the Hamiltonian becomes 


-iecAn-R 


n' = u-^nu = 


U = e 


{P - qAn - 2ecApf {p-qwA 


2M 


+ 


2mr 


+ H(|Pi-P2|) , 


where 


Qw — eiP 2 + ^ 2 Pi • 


The eigenfunctions of H and P,' are related via 

T' = Te*' 

The two cases of interest are 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


• ec ^ 0 , 

• q — 0 and P — 0 

Both lead to similar equations; so here we are going to focus mainly on the first case. 
For Cc — 0, the eigenfunctions can be factorised; 



T(cR,p) = x(cR)^(p) , 


(24) 


n^ = {ER + Ep)m . 


(25) 

The equation for the CMF variable y renders Landau levels (e.g., [I2])- The equation for the relative motion 

acquires the form of 





fnr^lp^ TC/I IN m 

2mc 8 

4 (p)= 0 , 

(26) 

where 

^1^ 
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(27) 

The vanishing of the commutator [H , ^2] = 0 permits further factorisation: 



V'(P) = , 


(28) 

where $((p) is the eigenfunction of the relative angular momentum operator 



d)(V9) = e“^ , 


(29) 


4 <h = s<h , s = 0 ,±l,±2. 


(30) 

The ensuing equation for ( 

is: 



+ 4 - + V{p) - e\ Cip) = 0 . 

2mr 2mrP p^ 2 8 

(31) 
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4 Various potentials 


C\C2 

In [12], the case of the Coulomb interaction V{p) = - between two charges was considered. Here we 

explore other potentials which permit exact solutions. 


4.1 Potential I 

We begin with 

Then the substitution 


r (p) ^ + .^ + fci p + ^2 p^ 

p p 


C,(p) = e-'’ -’’V*P.(P) 


(32) 

(33) 


furnishes the following equation for 


-pdp + (4rp^ + 2pp - 1 - dp + -p^ {ojlml + Sk 2 mr - 16r^) + p^ {2kimr - 4pr) + 


+p (4t{1 + - mr{sujc + 2Ep} - + 


— 4^^ + 4s^ 
4p 


(34) 


Ps = -[e + p(l + 20] Ps , 


where 


e = 2mreiei 


If the coefficients r, p and ^ are chosen to be 


T = - 


V = 


kiTUr 


2ki 


2 + ^ 
c mr 


2 r 

^ = \/ + 2'dm. 


(35) 

(36) 

(37) 

(38) 


the equation fl341) becomes: 


-pdl + (4rp^ + 2pp-l- 2^) dp + p (4r{l + O “ mr{suJc + 2Ep} - p^) + 


Ps = -[e + p{l + 20] Ps 


Our next goal is to cast the equation fl3^ into the form ([HD or, to be more exact, into the form 


'Jn^n + Jn ~ 


Ps = -UPs 


(39) 


(40) 


which is the special case of ((HD. We remind that j“, Ot, are generators of an irreducible n—dimensional 
representation of the SL{2, C) group. 

To that end, we insert into fl^Up the operators in the differential form, and change the variable as cp —)■ p: 


-pdp + ( c p + -{n- 2a) + (3cp ) dp - nc p 


(3n 

Ps = - \ l^- ^ ] CPs 


(41) 
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Equating of the coefficients in 


to those in fHTl) gives us: 

n 
2 

AkiTTlr 


ol — 1 + 2(^ + 


/? = - = - 

c {mlul + Sk2mrf^'^ 


c = 2^/T = ^mlul + 8 /C 2 "ir 
The energy of the relative motion for a fixed n is: 


(42) 

(43) 

(44) 


Ep = 


4r((^ + 1 ) — sUcfUr + c^n — rf 

2mr 


U)r 


Sk 


1 H-+ 1 + \/s^ + 2'dmr) — s — 


Akl 


ujlrrir 


^ I ^cmr 


We see that the infinite degeneracy of the energy for |s| > 0, which was the case in 
for nonzero "d and k 2 - 

By comparing RHS of fl39D and fITT]) . we find the condition 


e + 77(1 + 2 ^) = c ( 12 - 


Inserting therein the expressions (jSTj), fl43|) . and 
2ki 


e + 


(An 

T 

for rj, 13 and c from, we arrive at 

AkiTTlr 


(45) 

, does not take place 


(46) 


^ 


(1 + 2\/s^ + 2'i3mr) = ^mlul + 8k2mr ( — 


n 


2 {mlul + Sk2mr) 


3/4 


(47) 


This equation serves as a constraint upon the permissible values of the parameters entering our problem. 
(These include the magnetic field, along with the constants in the expression fl5^ for the potential.) Permis¬ 
sible is a set of values of these parameters, for which the equation (140]) . or its equivalent fITT]) . has a solution. 
The values of all these parameters, except one, can be chosen arbitrarily. The remaining parameter will then 
assume a discrete set of values, to satisfy the quantisation condition (147)) . With the potential fixed, this may 
be interpreted as ‘quantisation’ of the magnetic field. This interpretation will work also for the next potential 
considered. 

In case ki, /c 2 and d are equal to zero,0 the above equation reduces to 


7^2 = 


UcTTlr 


(48) 


which coincides with the equation (25) in 

At this point, two questions emerge. First, what if the expression under the square root in (145)) is negative? 
Second, how to describe the fall on the centre due to the presence of the inverse-square term in the potential? 
The two issues turn out to be closely connected, and can thus be resolved simultaneously. 

For p —)■ 0 , the equation for C,{p) assumes the form of 


ap + Ap-^)C(p) = o . 


(49) 


^ Be mindful that vanishing of fci entails the vanishing of /3 and rj. 
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where 


A = + 2'dmr . (50) 

As is explained in greater detail in [15] , no fall towards the centre takes place for A ^ 0. This helps us to 
observe that the condition 

+ 2^mr ^ 0 (51) 

is quite reasonable, given the physical content of the problem. 


4.2 Potential II 

Consider the potential of the form 

^(p) = ~ ^2 + ^4 + pQ . (52) 

A substitution helpful in this case is 

c.(p) = • (53) 


Insertion thereof into the formula flTI]) entails the following equation for ps: 

—pdp + (Srp'^ + 4pp^ — 1 — 2,^) dp + ^p^ {Skinrir — 64pr) + ^p® (Sk^mr — 64r^) + 
+ ^p^ ( 8 k 2 mr + ojlml - 16 {rf - 2(^ + 2)r) j + p (4p(^ + 1) - rnr{ 2 E + suic)) + ^^^— 2 — 


wherefrom we hud: 


T = 


\/2kQmr 

4 




k4mr y'2mrk4 

St 4v^ 

^ = \/s2 + 2Pmr 


•) 


(54) 

(55) 

(56) 

(57) 


Then the equation takes form 


-pdp + (Srp^ + 4pp^ - 1 - 2 ^) dp + ^p^ (sk 2 mr + ulml - 16 - 2(^ + 2)r) j + 


+p (4p(^ + 1) — mr{2E + sWc)) 


Ps = 0 , 


(58) 


We wish to cast the above equation fl55D into the standard form (1401) . To accomplish this, we substitute the 
operators in the differential form (the variable in fl40|) now being called p), and change the variable as cp —)■ p: 


-Pdj + 


c^p^ + c^p - a + ^ j - 


(?np 



(59) 
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After we change the variable one more time, p ^ the equation gets simplihed to 


-pd^p + (2c^p^ + 2j3cp^ — 2 a + n + 1 ) dp — Ac^np^ + 2c{n(3 — 2 u)p 


Ps = 0 . 


It acquires the needed form if we choose the parameters as: 

ot = I + i + — , 

R - JL - h ( 
c = 2 ^/t = { 2 kQmrY^^ ■ 

The energy of the relative motion will be: 

c(n/3 — 2 u) 2 p 


Ep = 


nir 


+ -(1 + 0-^ , 

nir 2 


while the quantisation condition will read as 

1 


= 4?7^ — 4c^n — 8t{^ + 2) — 2k2mr , 


— a/ 2kQmr (^An + 2^/+ 2'dmr + 2j — 2^2w-r 


2k 
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5 Potential III 

Finally, let us explore the option 


V{p) — — + ^ + — -\ - k2p 

p^ P'^ p^ p 


This time the needed substitution and the equation for ps are: 

C(p) = e~^~^^p^ps{p) , 

-p^dp + ( 2 pp^ - (2^ + l)p - 2r) dp + - 8 k 2 mr) - p^ ( 2 Emr + + sUcTrir) + 


+p(2,« + 2l,ra,. + „) + - 2{r + r 2_Un^ 


P 


P- 


Ps = - 2pr - - 2l2mr) Ps 


We compare its coefficients with 

-J-J+-2aJ- + /?J° + 27J^ 

By substituting operators in the differential form, we get 

+ ( 27 P^ + RP + n- 2 ) - 2a) dp - 2 'pip 


p, = -vp. 


Pn 

p, = — + n|p. 


(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 


(67) 

( 68 ) 

(69) 

(70) 


(71) 



From where we find: 


a 


-f = r] = 


= T = \/2Umr 
llTTlrT 


IsTUr + T{n + 1) 


^ 1 Unir 

«=2 + V 

21^171^ 

p = - n 


(72) 

(73) 

(74) 

(75) 


The energy of the relative motion is 


1 / r7‘ 

Eq = -I SUJr - 

^ 2 I 2m, 


SOJr. 


hy/2Umr 


2(l^mr + \/2/4m,{n + 1}) 


The hrst constraint on the parameters of the potential has the form of 

2/2m, — 2r]T — ^^ + s^ = i/ + n(l — ^ 


or, after we snbstitnte the constants: 

4/i/4m^ 


2 / 2 ^, — 


/ 3 m, + \/2l4r?v(n + 1 ) v 2 \/2UmrJ 


1 hrur Y 


I 2 I / 1 I n 

+ s = u + n\ IH- , . + 


y/2Umr 2 


(76) 


(77) 


In this case, the qnantisation of parameters of the potential cannot be rednced to the qnantisation of the 
magnetic held. Instead, we have one more constraint: 


w^m, = 8^2 . 


(78) 


6 The neutral system 

Let ns now discuss the case of ei = —62 = e and show that under certain circumstances it also renders 
exact solutions. Although in the general case the variables in the unitary-transformed Hamiltonian 


^ ^ (P-eB xpf ^ (p - ^ _ ^^1) 


(79) 


2M 2m, 

cannot be separatedjl we nevertheless can consider i?-independent solutions corresponding to a composite 
particle at rest: 

vl/'(i^,p) = ^(p) . (80) 

The equation for the relative motion takes the form 


p 

2 m, 


1 . ^ mrnlp 

2-V^ 


uJa^. + —^ + V{\p\)-Ep 


i){p) = 0 , 


(81) 


^ It is worth noting that in some situations, where the variables in the Hamiltonian cannot be separated, the wave function 
can nevertheless be factorised. 
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where 


(82) 


nn = 


eB 

2 mr 


and 


Uq = 


eB |/i 2 - yWlI 


nir 


(83) 


Further steps are analagous to the case of Cc = 0. The factorisation (1281) leads to the following equation 
for C(p): 


^ -d^ - dp + ^ - + 




2mr ^ 2mrp 2 


+ V{p) - Ep 


C(P) = 0 


(84) 


which is almost identical to fl3T]) . Below we present an inventory of changes in the solutions for the three 
afore-considered potentials. In what follows, Uq always replaces Uc- 

6.1 Potential I 

For the potential of the form ([32]), we again use the substitute (l33D . changing this time fl36D to 

r = Sk 2 mr . (85) 


The equations (15^ - (1441) remain unchanged. The energy of the relative motion and the quantisation condition 
will look as 


Ep = 


4:t{^ + 1) — sVtqnir + c^n — rf 1 


2 mr 


\ AVtl + —(n + 1 + \/s‘^ + 2'dnnr) - sUq - -j 
V nir I 


Ak\Uq 


Ani + ^ ) rur 

Q nir J ' 


e-1- 


2 ki 


4^2 + 

Q mr 


(1 2 \/s^ + 2 'dmr) = + 8 k 2 mr [ ^ ^ 


Akifrir 


(4m2f22 _|_ Sk2mryB 


( 86 ) 

(87) 


6.2 Potential II 

In this case, we still employ the equations 




rUr rUr 

The quantisation condition assumes the form of 

rrirk"^ 


and arrive at the same expresssion for the energy: 

( 88 ) 


^2 2 ^ 

^ " 2 kf, 


— \/2k^m^ {Au -|- 2^5^ -|- 2'dmr -|- 2 j — 2k2mr 


(89) 


6.3 Potential III 

In this situation, the equations 
the hrst constraint: 


m 1 / V 

Ep = \ SUq + 


(1751) stay unchanged, and so do the expression for the energy and 

hr 


2 mr 


SUq + 


2 {hmr + r{n 1}) 


(90) 
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2 


(91) 


2l2'mr — 


4Zi^4m: 


1 Isnir 

2 


l^rrir + \/2l^m^{n + 1 ) 

The second constrain, however, is now subject to change: 


+ s = ly + n I 1 + 


l^rUr n 
\/2l^rn^ 2 


VL^^rUr = 2k2 


(92) 


7 Conclusions 

In the presented paper, we extended the class of quasi-exactly solvable problems. We considered a problem 
of two charged particles placed in an external magnetic held and interacting with each other through a 
(generally, non-Coulomb) potential. The hrst of addresseed potentials was a generalizstion of the potential 
studied in [12]. All three considered Hamiltonians resemble those explored in n II They share a common 
feature: if all of the Hamiltonian’s parameters, except one, are hxed, then the exact solutions will exist for 
an inhnite discrete set of values of the remaining parameter. For the hrst two Hamiltonians, this observation 
may be written down as ‘quantisation’ of the magnetic held, with the potential’s parameters hxed. In the 
third Hamiltonian, however, the sole permissible value of the magnetic held is unambiguously linked to one 
of the hxed parameters; so this time it is one of the parameters of the potential that becomes subject to 
‘quantisation’. 

In all three cases, the energy depends on the parameters of the potential and magnetic held, as well as on 
the magnetic quantum number s and the dimension n of the representation of the generators Jn- Diagolalising 
a hnite n x n block of the Hamiltonian, one hnds n energy levels and the corresponding polynomial wave 
funtions. 
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